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Abstract 

The proton spectrum in the kinetic energy range 0. 1 to 200 GeV was measured by 
the Alpha Magnetic Spectrometer (AMS) during space shuttle flight STS-91 at an alti- 
tude of 380km. Above the geomagnetic cutoff the observed spectrum is parameterized 
by a power law. Below the geomagnetic cutoff a substantial second spectrum was ob- 
served concentrated at equatorial latitudes with a flux ~70m _2 sec _1 sr _1 . Most of these 
second spectrum protons follow a complicated trajectory and originate from a restricted 
geographic region. 
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Introduction 



Protons are the most abundant charged particles in space. The study of cosmic ray protons improves 
the understanding of the interstellar propagation and acceleration of cosmic rays. 

There are three distinct regions in space where protons have been studied by different means: 

• The altitudes of 30-40 km above the Earth's surface. This region has been studied with balloons 
for several decades. Balloon experiments have made important contributions to the understand- 
ing of the primary cosmic ray spectrum of protrons and the behavior of atmospheric secondary 
particles in the upper layer of the atmosphere. 

• The inner and outer radiation belts, which extend from altitudes of about 1000 km up to the 
boundary of the magneto sphere. Small size detectors on satellites have been sufficient to study 
the high intensities in the radiation belts. 

• A region intermediate between the top of the atmosphere and the inner radiation belt. The 
radiation levels are normally not very high, so satellite-based detectors used so far, i.e. before 
AMS, have not been sensitive enough to systematically study the proton spectrum in this region 
over a broad energy range. 

Reference Q includes some of the previous studies. The primary feature in the proton spectrum 
observed near Earth is a low energy drop off in the flux, known as the geomagnetic cutoff. This 
cutoff occurs at kinetic energies ranging from ~10MeV to -10 GeV depending on the latitude and 
longitude. Above cutoff, from -10 to -100 GeV, numerous measurements indicate the spectrum falls 
off according to a power law. 

The Alpha Magnetic Spectrometer (AMS) Q is a high energy physics experiment scheduled for 
installation on the International Space Station. In preparation for this long duration mission, AMS 
flew a precursor mission on board the space shuttle Discovery during flight STS-91 in June 1998. In 
this report we use the data collected during the flight to study the cosmic ray proton spectrum from 
kinetic energies of 0.1 to 200 GeV, taking advantage of the large acceptance, the accurate momentum 
resolution, the precise trajectory reconstruction and the good particle identification capabilities of 
AMS. 

The high statistics (~ 10 7 ) available allow the variation of the spectrum with position to be mea- 
sured both above and below the geomagnetic cutoff. Because the incident particle direction and mo- 
mentum were accurately measured in AMS, it is possible to investigate the origin of protons below 
cutoff by tracking them in the Earth's magnetic field. 

The AMS Detector 

The major elements of AMS as flown on STS-91 consisted of a permanent magnet, a tracker, time 
of flight hodoscopes, a Cerenkov counter and anticoincidence counters [3|. The permanent magnet 
had the shape of a cylindrical shell with inner diameter 1.1m, length 0.8 m and provided a central 
dipole field of 0.14Tesla across the magnet bore and an analysing power, BL 2 , of 0.14Tm 2 parallel 
to the magnet, or z-, axis. The six layers of double sided silicon tracker were arrayed transverse 
to the magnet axis. The outer layers were just outside the magnet cylinder. The tracker measured 
the trajectory of relativistic singly charged particles with an accuracy of 20 microns in the bending 
coordinate and 33 microns in the non-bending coordinate, as well as providing multiple measurements 
of the energy loss. The time of flight system had two planes at each end of the magnet, covering the 
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outer tracker layers. Together the four planes measured singly charged particle transit times with an 
accuracy of 120psec and also yielded multiple energy loss measurements. The Aerogel Cerenkov 
counter (n - 1.035) was used to make independent velocity measurements to separate low energy 
protons from pions and electrons. A layer of anticoincidence scintillation counters lined the inner 
surface of the magnet. Low energy particles were absorbed by thin carbon fiber shields. In flight the 
AMS positive z-axis pointed out of the shuttle payload bay. 

For this study the acceptance was restricted to events with an incident angle within 32° of the 
positive z-axis of AMS and data from two periods are included. In the first period the z-axis was 
pointing within 1° of the zenith. Events from this period are referred to as "downward" going. In the 
second period the z-axis pointing was within 1° of the nadir. Data from this period are referred to as 
"upward" going. The orbital inclination was 51.7° and the geodetic altitude during these two periods 
ranged from 350 to 390 km. Data taken while orbiting in or near the South Atlantic Anomaly were 
excluded. 

The response of the detector was simulated using the AMS detector simulation program, based on 
the GEANT package H. The effects of energy loss, multiple scattering, interactions, decays and the 
measured detector efficiency and resolution were included. 

The AMS detector was extensively calibrated at two accelerators: at GSI, Darmstadt, with helium 
and carbon beams at 600 incident angles and locations and 10 7 events, and at the CERN proton- 
synchrotron in the energy region of 2 to 14GeV, with 1200 incident angles and locations and 10 s 
events. This ensured that the performance of the detector and the analysis procedure were thoroughly 
understood. 

Analysis 

Reconstruction of the incident particle type, energy and direction started with a track finding proce- 
dure which included cluster finding, cluster coordinate transformation and pattern recognition. The 
track was then fit using two independent algorithms 0|6||. For a track to be accepted the fit was 
required to include at least 4 hits in the bending plane and at least 3 hits in the non-bending plane. 

The track was then extrapolated to each time of flight plane and matched with the nearest hit if it 
was within 60 mm. Matched hits were required in at least three of the four time of flight planes. The 
velocity, /3 = v/c, was then obtained using this time of flight information and the trajectory. For events 
which passed through the Cerenkov counter sensitive volume an independent velocity measurement, 
/3 C , was also determined. To obtain the magnitude of the particle charge, |Z|, a set of reference distri- 
butions of energy losses in both the time of flight and the tracker layers were derived from calibration 
measurements made at the CERN test beam interpolated via the Monte Carlo method. For each event 
these references were fit to the measured energy losses using a maximum likelihood method. The 
track parameters were then refit with the measured /3 and Z and the particle type determined from the 
resultant Z, /3, /3 C and rigidity, R = pc/\Z\e (GV). 

As protons and helium nuclei are the dominant components in cosmic rays, after selecting events 
with Z = +1 the proton sample has only minor backgrounds which consist of charged pions and 
deuterons. The estimated fraction of charged pions, which are produced in the top part of AMS, 
with energy below 0.5 GeV is 1 %. Above this energy the fraction decreases rapidly with increasing 
energy. The deuteron abundance in cosmic rays above the geomagnetic cutoff is about 2%. To 
remove low energy charged pions and deuterons the measured mass was required to be within 3 
standard deviations of the proton mass. This rejected about 3 % of the events while reducing the 
background contamination to negligible levels over all energies. 
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To determine the differential proton fluxes from the measured counting rates requires the ac- 
ceptance to be known as a function of the proton momentum and direction. Protons with different 
momenta and directions were generated via the Monte Carlo method, passed through the AMS de- 
tector simulation program and accepted if the trigger and reconstruction requirements were satisfied 
as for the data. The acceptance was found to be 0.15 m 2 sr on average, varying from 0.3 to 0.03 m 2 sr 
with incident angle and location and only weakly momentum dependent. These acceptances were 
then corrected following an analysis of unbiased trigger events. The corrections to the central value 
are shown in Table [T] together with their contribution to the total systematic error of 5 %. 



Correction 


Amount 


Uncertainty 


Trisser 

DO 






4-Fold Coincidence 


-3 


1.5 


Time of Flight Pattern 


-4 


2 


Tracker Hits 


-2 


1 


Anticoincidence 





1 


Analysis: 






Track and Velocity Fit 


-2 


1.5 


Particle Interactions 


+ 1 


1.5 


Proton Selection 


-2 


2 


Monte Carlo Statistics 





2 


Differential Acceptance Binning 





2 


Total 


-12 
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Table 1 : Acceptance corrections and their systematic uncertainties, in percent 

To obtain the incident differential spectrum from the measured spectrum, the effect of the detector 
resolution was unfolded using resolution functions obtained from the simulation. These functions 
were checked at several energy points by test beam measurements. The data were unfolded using a 
method based on Bayes' theorem ffZlOO, which used an iterative procedure (and not a "regularized 
unfolding") to overcome instability of the matrix inversion due to negative terms. Fig.[l]compares the 
differential proton spectrum before and after unfolding in the geomagnetic equatorial region, defined 
below. 

Results and Interpretation 

The differential spectra in terms of kinetic energy for downward and upward going protons integrated 
over incident angles within 32° of the AMS z-axis, which was within 1° of the zenith or nadir, are 
presented in Fig. Eland Tables[2j-|4l The results have been separated according to the absolute value of 
the corrected geomagnetic latitude J9), M (radians), at which they were observed. Figs. [2^., b and c 
clearly show the effect of the geomagnetic cutoff and the decrease in this cutoff with increasing M . 
The spectra above and below cutoff differ. The spectrum above cutoff is refered to as the "primary" 
spectrum and below cutoff as the "second" spectrum. 
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I. Properties of the Primary Spectrum 



The primary proton spectrum may be parameterized by a power law in rigidity, <J> x R r . Fitting [ 8 1 
the measured spectrum over the rigidity range 10 < R < 200 GV, i.e. well above cutoff, yields: 



7 = 2.79 ± 0.012 (fit) ± 0.019 (sys), 



4> = 16.9 ± 0.2 (fit) ± 1.3 (sys) ± 1.5 (7) 



GV 2 - 79 



m 2 sec srMV 

The systematic uncertainty in 7 was estimated from the uncertainty in the acceptance (0.006), the 
dependence of the resolution function on the particle direction and track length within one sigma 
(0.015), variation of the tracker bending coordinate resolution by ± 4 microns (0.005) and variation of 
the selection criteria (0.010). The third uncertainty quoted for <J> reflects the systematic uncertainty 
in 7. 



II. Properties of the Second Spectrum 

As shown in Figs. Ek, b, c, a substantial second spectrum of downward going protons is observed 
for all but the highest geomagnetic latitudes. Figs. EH, e, f show that a substantial second spectrum 
of upward going protons is also observed for all geomagnetic latitudes. The upward and downward 
going protons of the second spectrum have the following unique properties: 

(i) At geomagnetic equatorial latitudes, M < 0.2, this spectrum extends from the lowest measured 
energy, 0.1 GeV, to ~6 GeV with a flux -VOm^sec^sr" 1 . 

(ii) As seen in Figs. EH, d, the second spectrum has a distinct structure near the geomagnetic equa- 
tor: a change in geomagnetic latitude from to 0.3 causes the proton flux to drop by a factor of 
2 to 3 depending on the energy. 

(iii) Over the much wider interval 0.3 < 0m < 0.8, the flux is nearly constant. 

(iv) In the range < 0m < 0.8, detailed comparison in different latitude bands (Fig. |3j) indicates 
that the upward and downward fluxes are nearly identical, agreeing within 1 %. 

(v) At polar latitudes, M > 1.0, the downward second spectrum (Fig.Et) is gradually obscured by 
the primary spectrum, whereas the second spectrum of upward going protons (Fig.|2f) is clearly 
observed. 

To understand the origin of the second spectrum, we traced fTOl back 10 5 protons from their 
measured incident angle, location and momentum, through the geomagnetic field lITTIl for 10 sec 
flight time or until they impinged on the top of the atmosphere at an altitude of 40 km, which was 
taken to be the point of origin. All second spectrum protons were found to originate in the atmosphere, 
except for few percent of the total detected near the South Atlantic Anomaly (SAA). These had closed 
trajectories and hence may have been circulating for a very long time and it is obviously difficult to 
trace back to thier origin. This type of trajectory was only observed near the SAA, clearly influenced 
by the inner radiation belt. To avoid confusion data taken in the SAA region were excluded though 
the rest of the protons detected near the SAA had characteristics as the rest of the sample. Defining 
the flight time as the interval between production and detection, Fig. |4] shows the distribution of 
momentum versus flight time of the remaining protons. 
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As seen in Fig. H the trajectory tracing shows that about 30 % of the detected protons flew for 
less than 0.3 sec before detection. The origin of these "short-lived" protons is distributed uniformly 
around the globe, see Fig. |5^, the apparent structure reflecting the orbits of the space shuttle. In 
contrast, Fig. \5]p shows that the remaining 70 % of protons with flight times greater than 0.3 sec, 
classified as "long-lived", originate from a geographically restricted zone. Fig. |6] shows the strongly 
peaked distribution of the point of origin of these long-lived protons in geomagnetic coordinates. 
Though data is presented only for protons detected at M < 0.3, these general features hold true up 
to ®m ~ 0.7. Fig.|7]shows the distribution of the number of geomagnetic equator crossings for long- 
lived and short-lived protons. About 15 % of all the second spectrum protons were detected on their 
first bounce over the geomagnetic equator. 

The measurements by AMS in near Earth orbit (at 380km from the Earth's surface), between the 
atmosphere and the radiation belt, show that the particles in this region follow a complicated path in 
the Earth's magnetic field. This behavior is different from that extrapolated from satellite observations 
in the radiation belts, where the protons bounce across the equator for a much longer time. It is also 
different from that extrapolated from balloon observations in the upper layer of the atmosphere, where 
the protons typically cross the equator once. A striking feature of the second spectrum is that most of 
the protons originate from a restricted geographic region. 
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(20.6+ 1.2) xlO" 3 


(34.6+1.5)xl0" 3 

C91 9 + 1 91x1 O -3 

(116.+6.8)xl0- 4 
(57.2+4.7) xlO" 4 


(24.4+1. 1)xl0" 3 
C1 55 +Q 31x10~ 4 

(84.9+6.5)xl0" 4 
(40.0+3.8) xlO" 4 


(21.1 + 1.2)xl0- 3 

H91 +Q ^1x10~ 4 

(61.5+5.6)xl0" 4 
(26.9+3.4)xl0" 4 


(20.2+ 1.2) xlO" 3 
n 1 ^ +q oixi o~ 4 

(50.0+6.4) xlO" 4 
(24.2+4.2) xlO" 4 


1.15- 1.54 
1.54- 2.02 
9 09 - 9 69 
2.62- 3.38 


(116.+6.9)xl0- 4 
f66 Q+4 91x10~ 4 

(28.6+ 1.9) xlO" 4 
(110.+9.6)xl0- 5 


(28.6+3.3)xl0" 4 
C1 9 9 + 9 11x1 O -4 

^lz,,z,±z,,l y^lW 

( 8.2+1.8)xl0" 4 
( 3.6+1. 1)xl0 -4 


(17.7+2.5)xl0" 4 
( 8 5 + 9 61x10~ 4 

( 5.0+1.3)xl0" 4 
(30.0+8.6)xl0" 5 


(12.7+2.9)xl0" 4 
( 6 Q + 1 41x10~ 4 

(37.3+3.3)xl0" 5 
(204.+7.4)xl0" 6 


( 8.5±1.4)xl0- 4 
C 5 7 + 1 OlxlO -4 

(34.2+ 1.5) xlO" 5 
(29.0+ 1.4) xlO" 5 


3.38- 4.31 
4.31- 5.45 

5 45 6 86 
J.t-J — u.ou 

6.86- 8.60 


(44.3+7.9)xl0" 5 
nc 7+5 1 \ X 1 o~ 5 

( 6.1+2.2)xl0" 5 
(23.7+2.1)xl0" 5 


(20.3+6.0)xl0" 5 

n 3 4 + 4 81x1 0~ 5 

(105.+8.7)xi0" 6 
(53.8+2.7)xl0" 5 


(23.2+3.6)xl0" 5 

H7 6+3 91x10~ 5 

(31.9+2.3)xl0" 5 
(19.5 + 1.5)xl0" 4 


(25.0+1.3)xl0" 5 
f5S 5 + 5 Q1x10~ 5 

(32.1 +3.0) xlO" 4 
(96.2+6.4) xlO" 4 


(10.7 + 1. lJxlO" 4 

f69 Q + 6 41x1 O -4 
(18.4±1.4)xl0" 3 
(23.3 + 1.2)xl0" 3 


8.60- 10.73 
10.73- 13.34 
13 "KA.- 16 55 

16.55- 20.48 


(138.+6.8)xl0" 5 

C4Q 5 + 1 8">x10~ 4 

(65.7+2.1)xl0" 4 
(45.7 + 1.7)xl0" 4 


(28.6±1.7)xl0" 4 

f60 Q + 9 41x1 0~ 4 

(63.4+1.8)xl0" 4 
(45.5+1.7)xl0" 4 


(58.5+3.3)xl0" 4 

f85 7 + 3 11x10~ 4 
(72.1+2.1)xl0" 4 
(44.4+ 1.5) xlO" 4 


(128.+5.4)xl0" 4 
C1 1 5 +9 S1x10~ 4 

(75.6+2.5)xl0" 4 
(45.2+1.8)xl0" 4 


(193.+5.1)xl0" 4 

H9S + 1 ? 71X10 -4 

(75.6+2.7) xlO" 4 
(43.3 + 1.2)xl0" 4 


20.48- 25.29 
25.29- 31.20 
"X\ 90- 38 43 

38.43- 47.30 


(27.7 + 1.0)xl0" 4 
n 55 +5 Q1x10~ 5 

(90.5+4.1)xl0" 5 
(51.4+2.2)xl0" 5 


(25.5+1.0)xl0" 4 

C1 47 +7 1 1x1 0~ 5 
(79.2+4.7) xlO" 5 
(48.9+3.0)xl0" 5 


(255.+9.8)xl0" 5 

C144 +6 81x10~ 5 

(80.5 +4.5) xlO" 5 
(48.2+2.5) xlO" 5 


(248.+9.6)xl0" 5 

C149 +6 71x10~ 5 

(80.0+4.3) xlO" 5 
(48.2+3.0)xl0" 5 


(24.0±1.0)xl0" 4 

n +5 61x1 0~ 5 
(77.1+4.3)xl0" 5 
(47.1+2.7)xl0" 5 


47.30- 58.16 
58 16- 71 48 
71.48- 87.79 
87.79-107.78 


no 0+1 7ixio~ 5 
(164.+8.8)xl0" 6 
(86.1+3.9)xl0" 6 
(49.4+2.9) xlO" 6 


f98 6+9 01x10~ 5 
(15.4+1.2)xl0" 5 
(79.6+4.7) xlO" 6 
(45.0+4.6) xlO" 6 


C98 7 + 1 81x10~ 5 
(15.6+1.2)xl0" 5 
(81.5+6.4)xl0" 6 
(46.6+4.8) xlO" 6 


C9S 4+1 81x10~ 5 
(154.+8.8)xl0" 6 
(80.2+5.9) xlO" 6 
(45.8+2.8)xl0" 6 


C97 7 + 1 81x10~ 5 
(149.+9.9)xl0" 6 
(76.7+5.1)xl0" 6 
(43.4+2.6) xlO" 6 


107.78-132.27 
132.27-162.29 
162.29-199.06 


(28.6+3.1)xl0" 6 
(16.2+1. 8)xl0" 6 
(97.2+5. 1)xl0" 7 


(25.7+6.1)xl0" 6 
(14.3+7.0)xl0" 6 
(84.8+6.7)xl0" 7 


(26.9+7.3) xlO" 6 
(15.2+5.2)xl0" 6 
( 9.1+2.3)xl0" 6 


(26.4+6.2) xlO" 6 
(14.9+7.9)xl0" 6 
( 8.9±1.8)xl0" 6 


(24.8+4.6)xl0" 6 
(13.8+6.3)xl0" 6 
(82.1+6.2)xl0" 7 



Table 2: Differential downward proton flux spectra for lower latitudes. 
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Downward Proton Flux (m 2 sec sr MeV) 1 


Fr ■ 
^kin 

(GeV) 


Geomagnetic Latitude Range 


0.6 < @ M < 0.7 


0.7 < @ M < 0.8 


0.8 < @ M < 0.9 


0.9 < @m < 1.0 


1.0 < @M 


0.07- 0.10 
0.10- 0.15 
IS 99 

0.22- 0.31 


(12.2+3.5)xl0" 2 

( Q 7 + 1 ^x10~ 2 

(66.0+3.7)xl0" 3 
(44.4+ 1.6) xlO" 3 


(18.5±5.9)xl0" 2 
d 1 8+1 61x1 0~ 2 

(97.3±5.9)xl0 -3 
(44.2±2.0)xl0" 3 


(25.1 ±8.9) xlO" 2 
n q 1 +9 fiixio -2 

(144.±8.9)xl0" 3 
(92.4±6.9)xl0 -3 


( 4.3±1.3)xl0" 1 

C41 8 + 5 filxlO -2 

(33.6±3.3)xl0" 2 
(22.6±3.9)xl0 -2 


( 9.2±2.6)xl0" 1 
C q Q + i 9W10 -1 

(109.±6.7)xl0" 2 
(126.±5.3)xl0" 2 


0.31- 0.44 
0.44- 0.62 
o fo o 85 

0.85- 1.15 


(24.1 + 1.7)xl0" 3 
nos +s Rixio~ 4 

(47.8+6.7)xl0" 4 
(23.1 ±4.9) xlO" 4 


(23.8±1.3)xl0" 3 

C14 4+1 01x1 0~ 3 
(77.2±6.9)xl0 -4 
(60.9±6.5)xl0" 4 


(58.3±4.8)xl0" 3 
nfi 6+ 1 ? 5ixio~ 3 

(22.0±2.5)xl0" 3 
(34.9±5.8)xl0" 3 


(29.3±7.1)xl0" 2 
C a 7 + 1 n x 10 _l 

( 7.5±1.3)xl0" 1 
(85.3±7.5)xl0" 2 


(139.±4.1)xl0" 2 
n ^9 +4 sixi o~ 2 

(114.±4.2)xl0" 2 
(92.8±3.2)xl0" 2 


1.15- 1.54 
1.54- 2.02 
9 09 - 9 69 
2.62- 3.38 


(13.1±2.2)xl0" 4 
f 7 7 + 1 91x10~ 4 

(77.7±8.3)xl0" 5 
(49.1 ±5.9) xlO" 4 


(23.7±2.9)xl0" 4 

C44 8 + fi 71x1 O -4 

(43.1 ±5.8) xlO" 3 
(11.4±l.l)xl0- 2 


(15.4±2.4)xl0" 2 
C98 1 + 1 ? ^x10~ 2 

(30.9±1.8)xl0" 2 
(22.6±1.4)xl0" 2 


(71.7±4.5)xl0" 2 
^59 4+4 51x10~ 2 

(36.2±2.9)xl0" 2 
(24.8±2.1)xl0" 2 


(72.4±2.4)xl0" 2 
T51 1+1 41x10~ 2 

(37.0±l.l)xl0" 2 
(241.±6.4)xl0" 3 


3.38- 4.31 
4.31- 5.45 

5 45 86 
J.t-J — u.ou 

6.86- 8.60 


(27.9±2.9)xl0" 3 

f5fi 4 + 4 01x10~ 3 

(52.6±1.7)xl0" 3 
(35.6±1.2)xl0" 3 


(124.±4.6)xl0" 3 
f88 4 + 4 ^x1 0~ 3 

(55.6±3.2)xl0" 3 
(34.0±1.8)xl0" 3 


(15.4±l.l)xl0" 2 

fQ5 ^ + 5 QAxlO -3 

(59.3±3.5)xl0" 3 
(36.3±2.6)xl0" 3 


(16.2±l.l)xl0" 2 
nn - ! +7 7\xio~ 3 

(63.8±5.0)xl0" 3 
(39.0±2.8)xl0" 3 


(163.±3.1)xl0" 3 

H09 +9 QAx10~ 3 

(61.4±1.3)xl0" 3 
(390.±8.2)xl0" 4 


8.60- 10.73 
10.73- 13.34 

\"K "KA.- 1fi 55 
16.55- 20.48 


(212.±9.0)xl0" 4 

fiog _i_c o\„in-4 

(75.8±3.3)xl0" 4 
(41.7±1.5)xl0" 4 


(20.2±l.l)xl0" 3 

H91 +fi 41x10~ 4 

(69.0±3.8)xl0" 4 
(40.5±2.1)xl0" 4 


(21.8±1.6)xl0~ 3 

/-190 +0 01x10~ 4 

(75.2±4.3)xl0" 4 
(40.2±3.0)xl0" 4 


(22.5±1.6)xl0" 3 

Hi 1+1 ^IxlO -3 

(78.0±5.7)xl0" 4 
(39.3±3.3)xl0" 4 


(223.±6.5)xl0 -4 
n ^6 +4 5ixi o -4 

(76.2±2.7)xl0" 4 
(39.6±1.3)xl0" 4 


20.48- 25.29 
25.29- 31.20 
^1 90- ^8 4^ 

38.43- 47.30 


(24.9±l.l)xl0" 4 
n ^4 +5 fiixio -5 

(75.1±4.0)xl0" 5 
(46.0±2.7)xl0" 5 


(22.7±1.3)xl0" 4 
n ^9 +s 7^x1 o~ 5 

(69.2±4.5)xl0" 5 
(44.7±2.8)xl0" 5 


(237.±8.0)xl0" 5 

H97 +f\ 41x10~ 5 

(61.5±5.7)xl0" 5 
(44.0±3.5)xl0" 5 


(23.8±2.0)xl0" 4 
n? ^+1 4ixio~ 4 

(78.0±8.8)xl0" 5 
(44.1 ±4.6) xlO -5 


(22.0±1.3)xl0" 4 
n 1 s +7 Qixi o~ 5 

(76.7±6.5)xl0" 5 
(47.7±3.7)xl0" 5 


47.30- 58.16 
58 16- 71 48 
71.48- 87.79 
87.79-107.78 


C97 0+1 81x10~ 5 
(14.6±1.2)xl0" 5 
(76.0±4.6)xl0" 6 
(43.5±5.8)xl0" 6 


C96 3+1 91x10~ 5 
(142.±9.9)xl0" 6 
(72.9±4.5)xl0" 6 
(41.5±3.0)xl0" 6 


C95 7 + 9 81x10~ 5 
(13.9±1.3)xl0" 5 
(71.7±6.4)xl0" 6 
(41.1±4.1)xl0" 6 


C97 0+9 61x10~ 5 
(14.3±1.5)xl0" 5 
(72.5±6.5)xl0 -6 
(40.3±6.3)xl0" 6 


C98 5 + 9 61x10~ 5 
(154.±9.8)xl0" 6 
(79.3±8.7)xl0" 6 
(44.8±7.9)xl0" 6 


107.78-132.27 
132.27-162.29 
162.29-199.06 


(25.2±4.5)xl0" 6 
(14.3±3.9)xl0" 6 
( 8.6±1.5)xl0" 6 


(23.9±4.4)xl0" 6 
(13.4±4.7)xl0" 6 
(80.6±4.3)xl0" 7 


(23.9±4.4)xl0" 6 
(13.6±6.5)xl0" 6 
( 8.2±1.3)xl0" 6 


( 2.3±1.2)xl0" 5 
(12.3±8.9)xl0" 6 
( 7.2±3.7)xl0" 6 


( 2.6±1.2)xl0" 5 
( 1.4±1.4)xl0 -5 
( 8.5±2.4)xl0" 6 



Table 3: Differential downward proton flux spectra for higher latitudes. 
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Upward Proton Flux (m 2 sec sr MeV) 



Ekin 

(GeV) 


Geomagnetic Latitude Range 


@ M < 0.2 


0.2 < © M < 0.3 


0.3 < @ M < 0.4 


0.4 < @ M < 0.5 


0.5 < @ M < 0.6 


0.07-0.10 
0.10-0.15 
0.15-0.22 
0.22-0.31 


(16.4+4.4)xl0" 2 
(10.9+1.4)xl0- 2 
(85.3+4.9)xl0" 3 
(84.8+3.8)xl0" 3 


(13.1+3.9)xl0" 2 
( 7.5 + 1.0)xl0" 2 
(48.1±3.5)xl0" 3 
(44.5+2.1)xl0" 3 


(12.6+3.5)xl0" 2 
(66.0+9.2) xlO" 3 
(42.7+2.8)xl0" 3 
(39.3+1.9)xl0" 3 


(14.7+4.1)xl0" 2 
( 7.7 + 1. 1)xl0" 2 
(42.2+2.8) xlO" 3 
(35.5 + 1.8)xlO" 3 


(15.8+4.7)xl0" 2 
( 8.7 + 1.2)xl0" 2 
(46.3+2.8)xl0" 3 
(34.6+1.5)xl0" 3 


0.31-0.44 
0.44-0.62 
0.62-0.85 
0.85-1.15 


(66.8+3.4)xl0" 3 
(48.4+2.7) xlO" 3 
(32.7+2.0)xl0" 3 
(20.2+l.l)xl0" 3 


(33.6+1.7)xl0" 3 
(20.3+1.2)xl0" 3 
(120.+8.6)xl0" 4 
(53.9+4.6)xl0" 4 


(25.4+1. 1)xl0" 3 
(136.+8.3)xl0" 4 
(76.4+5.6) xlO" 4 
(42.0+4.5) xlO" 4 


(21.4+1. 1)xl0" 3 
(124.+9.2)xl0" 4 
(61.9+6.1)xl0" 4 
(31.9+4.6)xl0" 4 


(21.0+1. 1)xl0" 3 
(97.6+8.1)xl0" 4 
(34.8+4.3)xl0" 4 
(17.9+3.3)xl0" 4 


1.15-1.54 
1.54-2.02 
2.02-2.62 
2.62-3.38 


(124.+7.1)xl0" 4 
(62.0+4.2) xlO" 4 
(25.9+ 1.8) xlO" 4 
(10.7 + 1.5)xl0- 4 


(34.8+4.4)xl0" 4 
(16.4+2.3)xl0" 4 
( 7.9+1.3)xl0" 4 
( 4.2+1.2)xl0" 4 


(14.7+1.8)xl0" 4 
(12.5+2.3)xl0" 4 
( 5.6+1. ^xlO" 4 
(29.9+8.7) xlO" 5 


(14.0+2.3)xl0" 4 
( 8.8 + 1.8)xl0" 4 
( 4.6+1.2)xl0" 4 
(38.3 + 10.)xl0" 5 


( 8.6+2.1)xl0" 4 
( 5.2+1.2)xl0" 4 
( 3.4+1. ^xlO" 4 
(25.9+9.6) xlO" 5 


3.38-4.31 
4.31-5.45 
5.45-6.86 


(29.7+5.7)xl0" 5 
(11.2+4.6)xl0- 5 
( 3.7+2.4)xl0" 5 


(15.6+8.3)xl0" 5 
( 6.4+4.2) xlO" 5 


(11.9+4.9)xl0- 5 
( 7.2+3.8)xl0" 5 


(13.4+5.7)xl0" 5 
( 6.4+3.3)xl0" 5 


( 9.4+3.7)xl0" 5 


Ejcin 

(GeV) 


Geomagnetic Latitude Range 


0.6 < @ M < 0.7 


0.7 < @ M < 0.8 


0.8 < @ M < 0.9 


0.9 < @ M < 1-0 




0.07-0.10 
0.10-0.15 
0.15-0.22 
0.22-0.31 


(23.1 ±6.8) xlO" 2 
(10.5 + 1.5)xl0- 2 
(58.1+3.8)xl0" 3 
(43.0+2. 1)xl0" 3 


(32.9+9.5) xlO" 2 
(15.4+2.3)xl0" 2 
(72.5+5.4)xl0" 3 
(44.8+3.4)xl0" 3 


( 3.8+1. ^xlO" 1 
(18.0+2.4)xl0" 2 
(91.9+6.2)xl0" 3 
(57.4+3.3)xl0" 3 


( 5.1 + 1.5)xlO" 1 
(25.5+4.1)xl0" 2 
(99.8+8.4)xl0" 3 
(54.0+4.9) xlO" 3 




0.31-0.44 
0.44-0.62 
0.62-0.85 
0.85-1.15 


(20.7 + 1. 1)xl0" 3 
(83.4+8.0)xl0" 4 
(27.3 +4.0) xlO" 4 
( 7.2+2.3)xl0" 4 


(21.7+1.9)xl0- 3 
(78.6+9.3)xl0" 4 
(18.4+3.2)xl0" 4 
( 4.9+1.9)xl0" 4 


(25.7+2.6) xlO" 3 
( 8.8+1.2)xl0" 3 
(17.9+4.8)xl0" 4 
( 7.4+4.2)xl0" 4 


(22.5 +2.9) xlO" 3 
( 8.8 + 1.7)xlO" 3 
(23.4+8.0) xlO" 4 
(12.6+5.1)xl0" 4 




1.15-1.54 
1.54-2.02 
2.02-2.62 
2.62-3.38 


( 4.0+1.3)xl0" 4 
( 3.0+1.4)xl0" 4 
( 1.7 + 1. 2) xlO" 4 
( 6.3+4.1)xl0" 5 


( 3.2+2.3)xl0" 4 
(11.6+7.2)xl0- 5 
( 7.7+7.4)xl0" 5 
( 4.8+3.8)xl0" 5 


( 2.5+1.5)xl0" 4 
( 1.3+1.2)xl0" 4 


( 9.1+4.0)xl0" 4 
(16.8+9.3)xl0" 5 




3.38-4.31 


( 2.0+1. 1)xl0" 5 











Table 4: Differential upward proton flux spectra. 
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e M < 0.2 
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Figure 1: The proton differential flux in the equatorial region. Open circles show the measured 
distribution, filled circles are the data after unfolding. 
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Figure 2: Flux spectra for a,b,c) downward and d,e,f) upward going protons seperated according to 
the geomagnetic latitude, © M , at which they were detected. 
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Figure 3: Comparison of upward and downward second spectrum proton at different geomagnetic 
latitudes. As seen, below cutoff, the upward and downward fluxes agree in the range < 0m < 0.8 
(see also Figs.Eb, e). 
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Figure 4: The interval between production and detection, or flight time, versus momentum from the 
back tracing of protons detected in the region © M < 0.3. 
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Origin Longitude (degrees) 



Figure 5: The geographical origin of a) short-lived and b) long-lived protons with p < 3 GeV/c. The 
dashed lines indicate the geomagnetic field countours at 380 km. 
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Figure 6: The point of origin of long-lived protons (@ M < 0.3, p < 3 GeV/c) in geomagnetic coordi- 
nates. 
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Figure 7: Number of times the back traced trajectory crosses the geomagnetic equator for a) short- 
lived and b) long-lived protons (0 M < 0.3, p < 3 GeV/c). 
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